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ABSTRACT 

Using  high  frequency  (12-22  GHz)  VLBA  observations  we  confirm  the  exis¬ 
tence  of  a  Faraday  rotation  measure  gradient  of  ~500  rad  m-2  mas-1  transverse 
to  the  jet  axis  in  the  quasar  3C  273.  The  gradient  is  seen  in  two  epochs  spaced 
roughly  six  months  apart.  This  stable  transverse  rotation  measure  gradient  is 
expected  if  a  helical  magnetic  field  wraps  around  the  jet.  The  overall  order  to 
the  magnetic  field  in  the  inner  projected  40  parsecs  is  consistent  with  a  helical 
field.  However,  we  find  an  unexpected  increase  in  fractional  polarization  along 
the  edges  of  the  source,  contrary  to  expectations.  This  high  fractional  polariza¬ 
tion  rules  out  internal  Faraday  rotation,  but  is  not  readily  explained  by  a  helical 
field.  After  correcting  for  the  rotation  measure,  the  intrinsic  magnetic  field  di¬ 
rection  in  the  jet  of  3C  273  changes  from  parallel  to  nearly  perpendicular  to  the 
projected  jet  motion  at  two  locations.  If  a  helical  magnetic  field  causes  the  ob¬ 
served  rotation  measure  gradient  then  the  synchrotron  emitting  electrons  must 
be  separate  from  the  helical  field  region.  The  presence  or  absence  of  transverse 
rotation  measure  gradients  in  other  sources  is  also  discussed. 

Subject  headings:  galaxies:  active  -  galaxies:  quasars:  individual  (3C  273)  - 
galaxies:  jets  -  radio  continuum:  galaxies  -  polarization 


1United  States  Naval  Observatory,  10391  W.  Naval  Observatory  Rd.,  Flagstaff,  AZ  86002-1149 

2National  Radio  Astronomy  Observatory,  P.O.  Box  O,  Socorro,  NM  87801 

3Kavli  Institute  of  Particle  Astrophysics  and  Cosmology,  Menlo  Park,  CA  94025 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

13  MAY  2005 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Faraday  Rotation  Measure  Gradients  From  A  Helical  Magnetic  Field  in 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

US  Naval  Observatory  3450  Massachusetts  Ave,  NW  Washington,  DC 
20392 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

11 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


-  2  - 


1.  Introduction 

How  jets  are  launched  by  massive  black  holes  remains  one  of  the  fundamental  unsolved 
issues  in  astrophysics.  Many  researchers  have  suggested  that  magnetic  fields  are  intimately 
involved  in  the  collimation  of  the  jet  and  could  determine  which  sources  have  prominent  jets 
and  which  do  not  (Meier,  Koide,  &  Uchida  2001;  Koide  et  al.  2002).  Blandford  (1993)  urged 
observers  to  search  for  Faraday  rotation  measure  (RM)  gradients  transverse  to  the  relativistic 
jet.  Such  a  gradient  is  expected  if  a  helical  magnetic  field  wraps  around  the  jet,  and  among 
other  effects,  creates  a  Faraday  screen.  For  a  source  with  a  jet  axis  in  the  plane  of  the  sky 
a  helical  field  wrapping  around  a  jet  produces  a  maximum  line-of-sight  components  along 
the  jet  boundary,  no  net  line-of-sight  magnetic  field  component  at  the  center.  Blandford 
realized  that  this  would  produce  a  gradient  in  the  observed  RM  across  the  jet.  Projection 
effects  for  inclined  jets  produce  an  offset  in  the  RM  (Asada  et  al.  2002)  while  preserving  an 
RM  gradient. 

Faraday  rotation  measure  gradients  in  3C  273  (Asada  et  al.  2002)  based  on  5-8  GHz 
VLBA  observations,  present  intriguing  evidence  for  one  source  in  which  we  may  be  able  to 
see  the  effects  of  helically  wound  fields,  and  measure  the  sense  of  rotation  of  an  accreting 
black  hole.  With  the  search  for  these  RM  gradients  in  mind,  we  have  reanalyzed  12-22  GHz 
VLBA  polarimetry  observations  for  two  epochs  on  3C  273  to  obtain  higher  angular  resolution 
transverse  to  the  jet  axis,  while  maintaining  good  sensitivity  to  the  faint  polarized  emission. 

Throughout  this  discussion,  we  assume  H0=71  km  s-1  Mpc-1  (Spergel  et  al.  2003),  VtM 
=  0.27,  and  0.73.  This  gives  a  scale  for  3C273  of  1  mas  =  2.52  pc. 


2.  Observations  and  Results 

The  observations  were  made  on  2000  January  27  (2000.07)  and  2000  August  11  (2000.61) 
with  the  ten  element  Very  Long  Baseline  Array  (VLBA)1.  Data  reduction  and  calibration  for 
the  2000.07  epoch  were  discussed  in  Zavala  &  Taylor  (2001).  We  re-examined  the  data  for 
2000.07  and  made  polarization  angle  maps  using  using  five  frequencies:  12.1,  12.6,  14.9,  15.4 
and  22.2  GHz.  The  22.2  GHz  maps  were  tapered  to  approximate  the  12.1  GHz  resolution, 
and  all  maps  were  convolved  with  a  restoring  beam  matched  to  the  12.1  GHz  beam.  The 
observational  setup  for  the  2000.07  observation  was  repeated  for  the  2000.61  epoch.  The 
same  data  reduction  and  calibration  procedures  of  2000.07  (Zavala  &  Taylor  2001)  were 
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applied  to  the  previously  unpublished  2000.61  data.  A  datacube  of  polarization  position 
angle  maps  ordered  in  A2  was  used  to  create  the  rotation  measure  maps.  Pixels  in  the  RM 
maps  were  blanked  if  the  polarization  position  angle  error  was  more  than  ±  10°  at  any 
frequency. 

Rotation  measure  images  for  the  two  epochs  are  presented  in  Figure  la  (2000.07)  and 
Figure  lb  (2000.61).  Between  the  two  epochs  the  region  immediately  southwest  of  the  Stokes 
I  peak  (i.e.  in  the  projected  direction  of  the  jet)  shows  a  change  in  RM.  At  2000.07  the  RM 
is  «  1800  rad  m-2  and  6  months  later  shows  a  south  to  north  gradient  with  an  RM  from 
—  1000  to  +1000  rad  m-2. 

Beyond  3  mas  from  the  Stokes  I  peak  (hereafter  this  region  is  referred  to  as  the  “jet”) 
the  RM  distribution  becomes  smoother  and  does  not  change  dramatically  between  the  two 
epochs.  The  mean  RM  for  the  jet  is  642  rad  m-2  for  the  2000.07  epoch,  and  652  rad  m-2 
for  2000.61.  An  RM  gradient  transverse  to  the  jet  is  apparent,  and  we  discuss  this  in  §3.1. 

In  Figure  2  we  show  the  RM  corrected  intrinsic  B  vectors  for  3C  273  under  the  assump¬ 
tion  that  the  source  is  optically  thin  from  12-22  GHz.  We  simply  took  the  Faraday  corrected 
E  vectors  and  rotated  them  by  90  degrees.  This  was  done  for  the  entire  jet,  and  for  the  core 
west  of  the  Stokes  I  peak.  This  does  not  take  relativistic  aberration  into  account  (Lyutikov, 
Pariev,  &  Blandford  2003;  Lyutikov,  Pariev  &  Gabuzda  2005)  which  could  change  the  orien¬ 
tation  of  the  B  vectors  by  up  to  20°.  The  jet  B  vectors  are  initially  oriented  roughly  parallel 
to  the  projected  jet  direction,  then  rotate  to  almost  a  N— S  orientation  6  mas  southwest  of 
the  Stokes  I  peak.  This  changing  B  vector  orientation  repeats  itself  as  we  move  beyond  6 
mas  from  the  Stokes  I  peak.  For  the  2000.61  epoch  we  detected  polarization  as  far  as  15  mas 
down  the  jet,  and  the  B  vectors  there  have  returned  to  a  direction  parallel  to  the  projected 
jet  axis. 


3.  The  Faraday  screen  in  3C  273 

In  Fig.  1  we  noted  the  change  in  the  core  RM  at  3  mas  southwest  of  the  Stokes  I  peak. 
What  causes  the  change  in  the  observed  RM?  1  mas  corresponds  to  a  projected  distance 
of  2.5  parsecs  at  the  redshift  of  3C273.  A  region  of  that  size  is  unlikely  to  change  over  a 
six  month  time  interval.  We  thus  reject  a  change  in  the  physical  conditions  of  a  foreground 
Faraday  screen  as  the  cause  of  the  changing  rotation  measure.  A  more  likely  explanation 
for  the  RM  changes  in  the  core  over  the  6  month  time  interval  is  a  relativistically  moving 
polarized  jet  component  seen  through  a  foreground  Faraday  screen.  Changes  in  the  observed 
rotation  measure  are  caused  by  jet  components  sampling  different  sight  lines  as  they  move 
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out  from  the  core  (Zavala  &  Taylor  2004).  The  structure  in  the  core  is  consistent  with  the 
small  scale  polarization  structure  reported  by  Attridge  et  al.  (2003). 

The  Faraday  screen  in  front  of  the  jet  appears  to  be  relatively  constant  in  time  as  the 
RM  does  not  significantly  change  over  the  6  month  interval  presented  here,  or  over  three 
years  (Zavala  &  Taylor  2001).  The  jet  Faraday  screen  is  also  spatially  uniform,  as  opposed 
to  the  sub-structure  seen  in  the  core  (Fig.  1).  The  more  uniform  Faraday  screen  begins  at 
a  projected  distance  of  approximately  10  parsecs  from  the  core  (Fig.  1).  This  may  be  the 
distance  from  the  central  engine  beyond  which  the  magnetic  field  in  the  Faraday  screen  has 
become  more  uniform  and  smooth. 


3.1.  Evidence  for  a  systematic  RM  Gradient 

We  now  examine  the  case  for  a  transverse  RM  gradient  as  suggested  by  Blandford 
(1993).  We  wish  to  take  a  slice  across  the  broadest  RM  distribution  in  the  jet  transverse  to 
the  jet  motion.  We  use  Kellermann  et  al.  (2004)  to  define  a  jet  direction  of  —112.7°  at  the 
location  of  the  broadest  RM  distribution  in  Figure  1.  In  Figure  1  we  show  the  locations  of 
the  slices  across  the  jet,  with  the  orientation  of  the  slice  perpendicular  to  the  line  connecting 
the  Stokes  I  peak  of  the  component  at  a  relative  RA=  —3  mas,  relative  5  =  —6.5  mas  and 
the  map  center.  Figure  3  shows  the  resulting  RM  of  the  slices  at  the  two  epochs:  2000.07 
by  the  dashed  line  and  2000.61  by  the  solid  line.  A  transverse  gradient  of  approximately 
500  rad  m-2  mas-1  is  visible  in  Figure  3,  with  almost  4  beamwidths  along  the  gradient. 
This  is  much  larger  than  the  gradient  of  35  rad  m-2  mas-1  across  almost  two  beamwidths 
obtained  with  the  coarser  resolution  of  Asada  et  al.  (2002).  This  clearly  demonstrates  the 
advantage  of  working  at  the  highest  frequencies  the  VLBA  is  capable  of,  consistent  with 
enough  sensitivity  to  polarized  emission.  The  gradient  expected  by  Blandford  (1993)  is 
clearly  present,  and  is  suggestive  of  a  helical  field  wrapping  around  the  jet. 

The  RM  is  predominantly  low  and  mostly  negative  along  the  southern  edge  of  the  jet 
in  3C273  (ignoring  a  few  small  and  suspect  patches  of  very  high  RM),  even  to  a  distance  of 
15  mas  (38  pc)  in  epoch  2000.61.  Along  the  northern  edge  of  the  source  the  RMs  are  larger 
and  predominantly  positive.  This  overall  RM  structure  is  consistent  with  an  ordered  helical 
field,  and  suggests  that  this  field  is  the  dominant  source  of  the  Faraday  screen  on  projected 
linear  scales  out  to  at  least  ^40  pc. 
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3.2.  An  Alternative  View  —  Evidence  for  localized  RM  enhancements 

In  Figure  4  we  present  the  12  GHz  fractional  polarization  along  the  slices  shown  in  Figure 
1.  Along  with  a  gradient  in  the  RM  we  also  see  a  gradient  in  the  fractional  polarization. 
This  indicates  that  the  polarization  is  more  ordered  as  we  proceed  from  south  to  north  along 
the  transverse  RM  gradient.  Consider  the  simplest  case  of  a  jet  in  the  plane  of  the  sky  with  a 
helical  magnetic  field  wrapped  around  the  jet.  The  fractional  polarization  should  be  highest 
at  the  jet  center,  where  the  RM  is  zero.  The  fractional  polarization  then  falls  off  as  we 
approach  either  edge  of  the  jet,  as  the  RM  depolarizes  the  underlying  emission  (Gardner  & 
Whiteoak  1966).  Severe  projection  effects  are  surely  present  in  3C  273,  yet  it  is  still  surprising 
to  see  that  the  highest  RM  in  the  jet  corresponds  to  high  fractional  polarization.  Could  the 
localized  RM  and  fractional  polarization  enhancements  be  the  site  of  an  interaction  of  the 
jet  with  the  ambient  thermal  gas  which  forms  the  Faraday  screen  (Bicknell  et  al.  2003)?  In 
this  case  the  rotation  of  the  B  vectors  in  Figure  2  to  approximately  perpendicular  to  the 
jet  motions  seen  by  Kellermann  et  al.  (2004)  could  be  a  compression  and  enhancement  of 
the  local  magnetic  field  due  to  a  collision  with  the  Faraday  rotating  thermal  gas.  If  such 
an  interaction  is  responsible  for  the  observed  RM  and  fractional  polarization  gradients,  then 
the  lack  of  free-free  absorption  of  the  jet  of  3C  273  may  set  interesting  limits  on  the  electron 
density  of  the  thermal  gas  (Zavala  &  Taylor  2004). 


3.3.  RM  distributions  in  other  sources 

Sources  with  jets  broad  enough  to  be  well  resolved  are  rare  (Pollack,  Taylor,  &  Zavala 
2003),  but  a  few  examples  do  exist.  The  jet  of  M87  is  resolved,  but  the  RM  data  cover  a 
relatively  small  area  (Zavala  &  Taylor  2002).  There  is  no  obvious  RM  gradient,  but  the 
sign  change  of  M87’s  rotation  measure  is  consistent  with  a  helical  magnetic  field.  Gabuzda, 
Murray,  &  Cronin  (2004)  report  rotation  measure  gradients  in  4  BL  Lac  sources  with  an 
angular  resolution  comparable  to  Asada  et  al.  (2002).  3C  147,  a  compact  steep  spectrum 
(CSS)  quasar,  has  a  relatively  broad  jet  with  an  RM  gradient  observed  at  8  GHz  (Zhang 
et  al.  2004).  The  quasars  B1611+343  and  B2251+158  are  well  resolved,  but  there  is  no 
evidence  of  a  transverse  RM  gradient  (Zavala  &  Taylor  2003).  BL  Lac  is  marginally  resolved 
and  no  gradient  is  apparent  (Zavala  &  Taylor  2003).  Another  quasar  with  a  resolved  jet, 
B0736+017,  fails  to  display  an  obvious  transverse  RM  gradient.  At  present  we  can  say  that 
transverse  RM  gradients  do  exist,  but  these  gradients  are  not  a  universal  feature. 

The  very  high  fractional  polarization  we  observe  rules  out  internal  Faraday  rotation  as 
a  cause  for  the  observed  RM.  Internal  Faraday  rotation  is  expected  to  cause  severe  depo¬ 
larization,  and  discontinuities  in  the  observed  position  angle  versus  A2  plots  (Burn  1966; 
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Tribble  1991).  If  a  helical  field  is  responsible  for  the  observed  RM  gradient,  this  implies  a 
segregation  of  the  synchrotron  emitting  electrons  from  the  helical  field  region. 


4.  Summary 

We  confirm  the  presence  of  a  transverse  rotation  measure  gradient  of  500  rad  m-2 
mas-1,  significantly  larger  than  the  gradient  observed  by  Asada  et  al.  (2002).  This  rotation 
measure  gradient,  as  first  suggested  by  Blandford  (1993),  is  expected  if  a  helical  magnetic 
field  wraps  around  the  relativistic  jet  of  3C273.  This  field  imposes  an  order  on  the  Faraday 
screen  which  we  see  as  predominantly  negative  RMs  along  the  south  edge  and  positive  RMs 
along  the  north  edge  in  the  inner  projected  40  parsecs  of  the  jet.  The  increase  in  fractional 
polarization  along  the  RM  gradient,  and  turning  of  the  B  field  vectors  from  parallel  to 
perpendicular  to  the  projected  jet  motion,  suggest  that  an  interaction  may  also  explain  the 
gradient,  although  in  this  case  the  overall  order  of  the  RM  distribution  has  to  be  explained 
as  a  coincidence,  or  some  other  mechanism  that  imposes  an  overall  order  to  the  surrounding 
magnetic  fields.  The  magnetic  field  responsible  for  the  Faraday  rotation  appears  to  have 
structure  on  sub-parsec  scales  within  3  milliarcseconds  (8  pc)  from  the  core  of  3C  273.  The 
relatively  smooth  RM  structure  of  the  jet  implies  a  well  ordered  (on  parsec  scales)  magnetic 
field  in  the  Faraday  screen. 

The  larger  RM  gradient  and  the  gradient  in  fractional  polarization  structure  in  3C  273 
revealed  by  these  observations  illustrate  the  advantages  of  observing  broad  jets  at  the  highest 
frequencies  at  which  the  VLBA  operates.  The  superb  resolution  and  polarization  sensitivity 
of  the  VLBA  from  12-22  GHz  enables  high  spatial  resolution  tests  of  theoretical  predictions 
which  are  key  to  understanding  the  physics  of  relativistic  jets.  A  convincing  test  of  the 
helical  magnetic  field  model  awaits  a  clear  predominance  of  RM  gradients  with  the  expected 
helicity  in  many  sources,  or  the  detection  of  a  polarized  counter-jet  with  gradients  in  opposed 
directions  (Asada  et  al.  2002).  If  a  helical  magnetic  field  is  responsible  for  the  observed  RM 
gradients  then  because  the  high  fractional  polarization  observed  rules  out  internal  Faraday 
rotation,  the  synchrotron  emitting  electrons  must  somehow  be  segregated  from  the  helical  B 
field  region. 
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Fig.  1. —  Rotation  Measure  image  for  3C  273  on  two  epochs.  Total  intensity  contours  at 
12.1  GHz  are  overlaid  with  contour  levels  starting  at  12  m.Jy/beam  and  6.6  m.Jy/beam  for 
epoch  2000.07  and  2000.61  respectively  and  increase  by  factors  of  2.  The  synthesized  beam 
is  drawn  in  the  bottom-left  corner  and  has  dimensions  of  2.0  x  0.9  mas  at  0°  for  both  epochs. 
The  line  indicates  the  location  of  the  slices  shown  in  Fig.  3. 
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Fig.  2. —  Rotation  Measure  corrected  magnetic  vectors  for  3C  273  for  our  two  epochs  overlaid 
on  total  intensity  contours.  Contour  intervals  are  the  same  as  for  Figures  la  and  b.  Lengths 
of  magnetic  vectors  are  proportional  to  the  polarized  intensity  with  scales  of  1  mas  =  100 
mjy  beam-1  for  both  epochs. 
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Fig.  3. —  Rotation  Measure  image  slices  perpendicular  to  the  jet  axis  for  both  epoch  2000.07 
(dashed)  and  2000.61  (solid).  The  region  of  the  slice  is  indicated  by  the  solid  line  in  Fig.  1. 
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